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Abstract. We propose a scheme to experimentally observe matter- wave interference 
in the time domain, specifically in the arrival-time or the timc-of-flight (TOF) 
distribution for atomic BEC Schrodinger-cat state represented by superposition of 
macroscopically separated wave packets in space. This is in contrast to interference in 
space at a fixed time observed in reported BEC experiments. We predict and quantify 
the quantum interference in the TOF distribution calculated from the modulus of 
the quantum probability current density (rather than the TOF distributions obtained 
from a purely classical or semi-classical treatment in many reported experiments). 
The interference and hence the coherence in the quantum TOF signal disappears in 
the large-mass limit. Our scheme has the potential to probe the validity of various 
other theoretical approaches (Phys. Rep. 338, 353 (2000)) of calculating the quantum 
arrival time distribution. 
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1. Introduction 

In recent times, laser cooling and trapping of atoms has become an area of active 
research [I]. The temperature of the cold atomic sample is one of its most important 
characteristics and several methods have been proposed and used for its determination. 
A well-known technique of measuring this temperature is the time-of-flight (TOF) 
method [2]. It is significant to mention that the first evidence for Bose Einstein 
condensate (BEC) was emerged from TOF measurements [3]. Most of the samples 
of cold atoms are initially prepared in magneto-optical traps and the atomic cloud is 
allowed for a thermal expansion after its release from the trap. These so-called time-of- 
flight (TOF) measurements are performed either by acquiring the absorption signal of 
the probe laser beam through the falling and expanding atomic cloud, or by measuring 
the fluorescence of the atoms excited by the resonant probe light. Most of the theoretical 
analyses of TOF measurements are as follows. To find the shape of the absorption TOF 
signal, one assumes to start with the initial Gaussian position and velocity distributions 
of atoms in the trapped sample. The initial probability distribution of finding an atom 
in the phase space volume element with coordinates (zq, vq) is given by 

@(zq, v )dz Q dv = 1 exp 



x 77; — \yTF? ex P ( ) dz dv (1) 



Here for simplicity, we consider the one-dimensional case. The Gaussian width a v of the 
velocity distribution is associated with the temperature T of the cloud by the relation 
al = kT/m, where m stands for the atomic mass and k is the Boltzmann constant. 
Using the Newton's equations for ballistic motion of a particle accelerated by the earth's 
gravitational field (in the vertical z-direction) , the velocity is obtained in terms of the 
time of flight as 

2 9 >' ' dt ~ t 2 
Substituting the above expression for vo from Eq.(|2]) in Eq.((T]), and then finally 
integrating over Zq, one can obtain the TOF distribution at an arbitrary distance z = H, 
given by 

1 {\gt\2al + aie)-Haie) 
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Figure Q] shows a typical TOF or arrival time distribution for cold sodium atoms. 
This kind of purely classical analyses are adopted in most of the discussions on TOF 
measurements where arrival time of atomic or sub-atomic particles is treated as an 
elementary well-defined, unique, and classical quantity. Also, the theoretical treatments 
of the TOF distribution that can be obtained using, for instance, the Green's function 
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Figure 1. (Color online) Classical TOF distribution @(t) given by Eq.Q is plotted 
for BEC of sodium atoms at a temperature T = 1 \iK. The detector is located at a 
distance z = H = — 1 cm and ctq = 1 fim. 



method [I] or any semiclassical method [5], however, are equivalent to the TOF 
distribution obtained by using Newton's equations for ballistic motion of particles [2]. 
The interpretations or theoretical analyses of the results of the various TOF experiments 
[SI IE] with molecular, atomic or sub-atomic particles where classical trajectories 
are inferred from Newtonian mechanics remain debatable, especially in the domain of 
small atomic masses and low temperatures where quantum mechanical effects should 
be significant and quantum TOF distribution can not be reproduced with classical or 
semi-classical analyses. 

Here we provide an example in the context of BEC matter-wave interference, where 
a quantum analysis for TOF is necessary. We propose a scheme for measuring the 
TOF distribution of a freely falling atomic BEC prepared in non-classical Schrodinger- 
cat states. The interference in the TOF distribution (or signal) can then be observed 
by taking a note or record of the particle counts over various tiny time-windows at 
a fixed detector location. This is different from the interference between two freely 
expanding BEC's observed [9] in space after a definite time of free fall of the condensates. 
Coherent splitting of BEC atoms with optically induced Bragg diffraction have been 
done experimentally pUl E]. The spatial coherence of a BEC is measured using 
interference technique by creating and recombining two spatially displaced, coherently 
diffracted copies of an original BEC 

As mentioned above, most of the experiments (particularly when matter-waves are 
associated with centre-of-mass motion or external motion of massive quantum particles) 
demonstrate matter-wave interference by showing the intensity variation at an extended 
region of detection space at a fixed time. In the present paper, in contrast, we predict 
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and quantify the matter-wave interference in the center-of-mass motion by calculating 
the time distribution of matter- wave arrival probability at some fixed spatial point. More 
specifically, we discuss here the BEC matter-wave interference in the TOF distribution 
or arrival time distribution since BEC as a source of coherent matter waves is already 
routinely demonstrated and thus may be an ideal candidate to show interference signal 
in the time domain (arrival-time distribution). We use here a particular quantum 
approach to calculate the TOF distribution and in our analysis we do not use at any 
point classical or semiclassical ingredients. We consider the free fall of matter-wave 
associated to quantum particles represented by an initial Schrodinger-cat state which 
is the linear superposition of two mesoscopically distinguishable Gaussian wave packets 
peaked around different heights viz., z = and z = —d along the vertical z-axis. 
Then after a certain height of free fall (evolution under the potential V = mgz) of the 
Schrodinger cat, we calculate the quantum TOF distribution at a given detector location 
z = H. During the free fall, the distinct superposed wave packets of the Schrodinger 
cat overlap or interfere in space, so it is natural to expect that they will also interfere 
in the time of fall showing an interference pattern in the quantum TOF distribution. 

We take this particular example of matter- wave interference in the discussion of quantum 
TOF distribution to pinpoint the necessity of a quantum analysis. So, the need for 
a quantum analysis of TOF distribution is not merely a conceptual but a practical 
issue, asking how to predict the TOF distribution using only classical and semi- 
classical ingredients in a purely quantum scenario like this (interference in the TOF 
distribution for quantum particles). Now, in spite of the emphasis of quantum theory 
on the observable concept, there is no commonly accepted recipe to incorporate time 
observables and their probability distributions in the quantum formalism, and there 
is considerable difficulty and debate over the issue of defining time (for example, 
tunneling time, decay time, arrival time) as an observable [121 Q2J HU H3 EE H3 
H~9l 1201 I2T1 1221 1231 EH [25]. Even for the simplest case of arrival time problem there 
is no unique way to calculate the probability distribution in the quantum formalism 
|17j . Despite this, many researchers have evidently not been discouraged from seeking 
an expression for the arrival time distribution (or the quantum TOF distribution) 
within a consistent theoretical framework. Several logically consistent schemes for 
the treatment of the arrival time distribution have been formulated, such as those 
based on axiomatic approaches [H], operator constructions [15], measurement based 
approaches [HU [19], trajectory models [21] and probability current density approach 
[HI UHl HH [201 [2H E2l E31 [2H [25]. We will use here probability current density 
approach to calculate the quantum TOF distribution which is logically consistent and 
also physically motivating. 

The main purpose of our paper is two- fold. First, our proposal to experimentally observe 
or quantitatively predict the matter-wave interference in the time domain, specifically 
in the TOF distribution is itself quite significant which has not been explored in the 
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current literature to the best of our knowledge. BEC as a source of coherent matter 
waves is already routinely demonstrated in spatial interference experiments, so BEC will 
also be an ideal candidate to show the interference in the TOF distribution. Second, 
we have just mentioned that there is an inherent nonuniqueness within the formalism 
of quantum mechanics for calculating the TOF or arrival time distribution. It remains 
an open question as to what extent these different quantum mechanical approaches 
[121 H31 HH H5l H6l H71 [HI [19l [2OI [2H [22I [23I [2II [25] for calculating the time distributions 
can be tested or empirically discriminated. Our proposal of measuring matter-wave 
interference in TOF distribution has the potential to empirically resolve ambiguities 
inherent in the theoretical formulations of the quantum TOF distribution. In this 
respect, it would be interesting if the prediction of BEC matter wave interference in 
the TOF distribution (calculated from different quantum approaches) be verified in 
actual experiments. 

2. Interference in the quantum time-of-flight distribution for Bose-Einstein 
Condensate 

We begin our analysis with the standard description of the flow of probability in quantum 
mechanics, which is governed by the continuity equation derived from the Schrodinger 
equation given by 



The quantity J(x, t) = j^(^V^/* — \I/*V\I/) defined as the probability current density 
corresponds to this flow of probability. In one dimension, the current density J(x, t) 
tells us the rate at which probability is flowing past the point x. So, interpreting the one 
dimensional continuity equation in terms of the flow of physical probability, the Born 
interpretation for the squared modulus of the wave function and its time derivative 
suggest that the arrival time distribution of the particles reaching a detector located 
at x = X can be calculated [IZl HH HH1 [201 I2T1 [221 1231 [2H [25] using the probability 
current density J(x,t). It should also be noted that J(x,t) can be negative, hence one 
needs to take the modulus sign in order to use the above definition. Our aim here 
is to derive an expression for the TOF distribution through the quantum probability 
current density for the atomic BEC representing the mesoscopic Schrodinger cat and 
showing the interference in TOF signal. Probability current density approach to the 
TOF distribution is also justified by the Bohmian model of quantum mechanics in 
terms of the causal trajectories of individual particles [21]. Although the Schrodinger 
probability current density is formally nonunique up to a total divergence term [20] , the 
current can be uniquely fixed if one calculate the current in the non-relativistic limit of a 
proper relativistic wave equation which provide appropriate spin-dependent corrections 
to it [231 EE]. We ignore this small spin-dependent contribution here in our present 
discussion, as the estimated magnitude of the spin-dependent current is roughly 10 5 to 
10 6 times smaller than the Schrodinger current. It was emphasized that the probability 
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current density approach not only provides an unambiguous definition of arrival time at 
the quantum mechanical level [HI [HI UHl [20l [2H [22l [231 EH [25] , but also adresses the 
issue of obtaining the proper classical limit of the TOF of massive quantum particles 
El ESI. 



Now, to keep the discussion concise, here we restrict ourselves to the case of one- 
dimensional motion, but our resulting conclusion does not depend on three-dimensional 
extension which will be straightforward as discussed later at the end of this Section. 
A magnetically trapped BEC as a source of coherent matter wave or atom laser, 
where a macroscopic number of atoms occupy the same ground state is now routinely 
available. After being released from the trapped Bose-Einstein condensates, matter 
waves fall freely due to the gravity. If the atomic beam is well collimated, we can 
use a one-dimensional Gross-Pit aevskii equation 0, [27] for the evolution of condensate 
wavefunction ty(z,t) with the gravitational potential, 

ih—V(z,t) = - ——^(z,t)+mgz^(z,t) + U \^(z,t)\ 2 V{z,t), (5) 

where \^(z, t)\ 2 provides the density profile of the BEC, m denotes the atomic mass, g 
the gravity acceleration, and Uq the inter-atomic interaction strength. In our present 
discussion we consider condensate of non-interacting bosons and we neglect [271 [28] the 
effects of inter-atomic interaction Uq on the freely falling condensate. In the BEC, the 
whole complex is described by one single wave function ty(z,t) (a macroscopic wave 
function of the condensate) exactly as in a single atom, and we can speak of coherent 
matter in the same way as of coherent light in the case of a laser. To show interference 
in the quantum TOF signal for the freely falling BEC, we consider the initial state of 
the BEC be prepared in a Schrodinger-cat state which is the coherent superposition of 
two mesoscopically distinguishable states in the configuration space: 

*(z, 0) = JV [ciVifo 0) + c^ 2 (z, 0)] , (6) 

where 

are the Gaussian wave packets centered around z = and z = —d, respectively, and <To is 
the initial position spread. A description of the initial ID wavefunction for two separated 
BEC's using the Gaussian form ([6]) has been made, for example, by the authors of [29J, 
where they consider all the non-interacting bosons are prepared to be condensed in the 
ground state of the harmonic trap [3] . For simplicity, we take c\ = c 2 — 1/ y/2 which 
implies that after coherent splitting of the original BEC, each component has equal 
number of atoms. Then the value of the normalization constant 



JV = + exp(-rf 2 /8a 2 ). (9) 
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As we have mentioned, the Schrodinger-cat state of matter was generated for a BEC 
represented by superposition of spatially separated states and the superposition was 
verified [H] by detecting the quantum mechanical interference (in space) between the 
localized wave packets separated by a mesoscopic distance. Under the experimental 
situtations (where the spatial coherence of the BEC was measured using interference 
technique by creating and recombining two spatially displaced, coherently diffracted 
copies of an original BEC) discussed by the authors of [TUl [TT] . the BEC wave function 
can be written as a linear superposition of spatially separated wave packets [29] which 
may be inferred as a true macroscopic Schrodinger cat. Now considering the free fall 
of the coherently splitted BEC under gravity, we calculate the time evolution of the 
Schrodinger-cat state (jSJ) according to equation (JSJ) with Uq — 0, we then obtain 



V(z,t) = ^=[il> 1 (z,t)+fo(z,t)] 



where 



ipx(z,t) = (2ns 2 ) 1 ^ 4 exp 
and 

exp 



u^t) = (2^y 



(*+iggy 

(z + d+ \gt 2 ) 
4:S t a 



exp 



2+3 



i(jr) [ Qtz+-gH 



exp 



gtz + -g 2 t 3 



(10) 



(12) 



with 



St = cr (l + iht/2mal) . (13) 

The expression for the Schrodinger probability current density corresponding to the time 
evolved state ^(z, t) (TTUI) is given by 

^ 2 



where 



[Jx(z,t) + J 2 (z,t) + J 3 (z,t) + J 3 *(z,t)} 



Ji(z,t) 



J 2 (z,t) 



hH 



i 



hH ,1 2N 



x |^iM)| 2 , 



gt 



where 



J 3 (z,t) + J 3 *(z,t) = 2 P 12 (z,t) (77 cos 5 - A sin S) 

hd 



(14) 

(15) 

(16) 
(17) 



A 



v 



4ma 2 ' 

hH 



(2z + d + gt 2 ) - gt, 



and 



8m 2 <To 2 a 2 

Pl2(z,t) = \^(z,t)\ \Mz,t)\ 



(19) 
(20) 
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with the time-dependent position spread given by 

a 2 = St s* t = al (1 + h 2 t 2 /4m 2 a$ . (21) 

Here the quantity s t is defined in equation (|T3|) . The oscillatory factor 5 in (fT7|) . 
responsible for the interference effect, is given by 

fit . o , , 2 „ a Ht (d 2 + dgt 2 + 2zd) 
8ma V 2 8m + f^) 

Now taking the modulus of the quantum probability current density (fT4j) . we obtain 
the quantum TOF distribution [13 HSl UHl [201 [2D 1221 1231 [2H at a detector location 
z = H for the spatially separated BEC Schrodinger cat falling freely under gravity given 
by 

Tl(t) = \J(z = H,t)\. (23) 

Exactly similar expression for quantum TOF distribution can be obtained for a three 
dimensional analysis of the problem by considering the form of the initial wave function 

jY 

y, z,0) = -^= [faix, y, z, 0) + if) 2 {x, y, z, 0)] , (24) 

^ 1(X '^' 0) = ( 2 Jp ieXP (~WJ ^ ("^) ^ ' (25) 

Mx, V, *, 0) = (^71 exp j exp j exp ) , (26) 

where if)i(x, y, z, 0) and ^(x, ?/, ^, 0) are now three dimensional Gaussian wave packets 
separated along the vertical z-axis, having peaks around the points (0, 0, 0) and (0, 0, —d) 
respectively. The value of the normalization constant JV remains the same as that in 
Eq.flHJ). One can then obtain the three dimensional Schrodinger time evolved wave 
function under the gravitational potential. The quantum TOF distribution can then be 
calculated again using the three dimensional quantum current. Interpreting again the 
three dimensional continuity equation (j3J) in terms of the flow of physical probability, 
one can define the quantum TOF distribution for the atoms crossing a surface element 
dS as |J.dS|. It is important to mention here that the quantum flux density |J.dS| have 
been identified with the "time distribution" of particles crossing the surface element dS 
by Daumer et al. [30], who were applying Bohm model to the scattering problem for a 
quantum particle in three dimension. Hence quantum TOF distribution for the atoms 
reaching at a finite surface plane S in three-dimension will be given by 



II(t) = | J j J.dS | = | J j J.n dS |, (27) 

where n is the unit vector normal to the surface. The quantum TOF distribution for 
the atoms reaching the XY plane (n = — z) after a certain height (z = H) of free fall is 
then given by 

Ili(t) = | / / J.ii dS | = | / J z (x,y,z = H,t)dxdy |, (28) 
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where J z (x, y, z = H, t) is the z-component of the three dimensional probability current 
density at a fixed height z = H. By evaluating the integral of Eq. (1281) . one can see 
that Ili(t) is exactly the same as the quantum TOF distribution obtained for one 
dimensional analysis (1231) . To understant more clearly the origin of this interference in 
TOF, let us consider the propagation (evolution) of individual wave packets ipi and ip2 
under the gravitational potential. Then one will have two distinct TOF distributions 
having separate mean arrival times. This is because the peaks of the component wave 
packets ipi and ip2 take different times to reach the detector at z = H, since they are 
spatially separated along the vertical z-axis, and the interference in TOF arises due to 
the superposition of these two wave packets. In this setup, the cross term (interfering 
term) in the quantum TOF distribution arises from the relative phase of the component 
wave packets (ipi and 1P2) along z-direction, as only z-componets of the component wave 
packets differ in the time evolution and continue to develope the relative phase, and 
this relative phase is not cancelled out when we perform the integration over XY-plane. 
The interference pattern in the quantum TOF distribution can be detected by using a 
probe laser, focused in the form of a sheet underneath the falling BEC atoms in the 
XY-plane at z = H. When the trapping forces are turned off, the BEC atoms will 
fall through the laser probe under the influence of gravity. It is then possible to detect 
the fluorescence from the atoms excited by the resonant probe light as they reach the 
detection sheet. The fluorescence can be measured as a function of time to determine the 
TOF distribution. For this setup, one can also consider a situation where the detection 
is made in the YZ-plane (n = — x) at a fixed x = X. In that case, quantum TOF 
distribution (say, U 2 (t)) can be obtained from the x-component of the three dimensional 
current (J x (x = X,y,z,t)) integrated over YZ-plane using Eq.(|27l). By evaluating that 
integral, one can see that there will be no interference in the quantum TOF distribution 
n 2 (t) under this situation. This is because the interference term in the quantum TOF 
distribution is wiped out when we perform the integration over the YZ-plane. 

Next, we consider another setup in three dimension where the superposed wave packets 
ipi(x, y, z, 0) and ^(x, y, z, 0) are separated along the horizontal X-axis, having peaks 
around the points (0, 0, 0) and (— d, 0, 0) respectively. This situation is analogous to the 
experimental setup of Ref. [9] where BEC interference was observed in space at a fixed 
time. For this geometry, one can again calculate the three dimensional Schrodinger time 
evolved wave function under the gravitational potential. The quantum TOF distribution 
( 1271) can then be calculated again using the three dimensional quantum current. For 
this setup, we again consider the detection of the particles at a surface plane (XY-plane 
with n = — z) at z — H. In this case, quantum TOF distribution (say, Il3(t)) can be 
obtained from the z-component of the three dimensional current integrated over XY- 
plane using Eq. fl27|) . By evaluating that integral, one can see that there will not be 
any interference at all in the quantum TOF distribution r±3(t). The interference term 
in the quantum TOF distribution is wiped out when we perform the integration over 
the XY-plane. This is because for this setup, the individual wave packets ip\ and ip2 
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are not separated along the vertical z-axis, so they will have the same TOF distribution 
with the same mean arrival time to reach the detection plane at z = H . Hence, we do 
not expect any interference in the quantum TOF distribution detected at the horizotal 
plane at z = H when we consider the superposition of the horizontally separated wave 
packets. The interference term in the quantum TOF distribution, in this case, is wiped 
out when we perform the integration over the XY-plane, even though one can observe 
the interference in space at a fixed time. For this setup, one can also consider a situation 
where the detection is made in the YZ-plane (n = — x) at a fixed x = X. In that case, 
quantum TOF distribution (say, n 4 (t)) can be obtained from the x-component of the 
three dimensional current integrated over YZ-plane using Eq. (127jl . For this situation, 
although we will see the presence of some interfering terms in the expression of quantum 
TOF distribution Il^t), the intensity to observe this interference will be very low, as 
only a small fraction of the condensate atoms will arrive at the detection YZ-plane at 
x = X due to free expansion (free particle motion) of the wave packets. One can also 
check that the quantum TOF distribution Tl^t) will have exactly the same expression 
as II(t) of fl23|) with g = (no gravity) and with z = H replaced by x = X. We will 
show numerically in the next section that gravity plays an important role in our setup 
to pull down the vertically separated superposed condensate toward the detection plane 
at z = H. 

Hence the only two situations (in our above discussion) where we see the presence of 
interference in three dimensions are the quantum TOF distribution Ili(t) and n 4 (t). 
Now, Ili(t) is exactly the same as n(t) and H^t) is also same as n(t) with g = and 
with z = H replaced by x = X. So the whole characteristic of the interference pattern 
in quantum TOF distribution hinges upon the form of n(t). In the next section we 
study numerically the parameter dependence of the quantum TOF distribution n(t) 
and the physical interplay between these parameters. 

3. Numerical Results and Discussions 

Quantum TOF distribution Tl(t) of the freely falling atomic BEC Schrodinger cat is 
plotted (Figj2j) at a detector location z = H = —1 cm with different values of wave 
packet separation d. We see clear signature of interference in the quantum TOF 
distribution arising due to the terms J^(z = H,t) and Jz*(z = H,t) of (JTTJ) in the 
expression for quantum probability current density fl23l) and ffl^l) . During free fall, the 
spatially separated wave packets of the BEC Schrodinger cat overlap or interfere in space 
and hence they also interfere in the time of fall showing an interference pattern in the 
quantum TOF distribution. The quantum TOF distribution H(t) may be visualized as 
a coherent pulse of BEC atoms. The interference pattern in the quantum TOF signal 
(Fig 12^) is very sharp for a typical set of parameter values, for example, H = — 1 cm, 
d = 50 /mi, o"o = 1 /^ m ; and the pattern disappears (FigEJ) when the separation between 
the BEC superposed wave packets is decreased to d = lfim for the above mentioned 
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Figure 2. (Color online) Quantum TOF distributions H(t) for the coherently splitted 
BEC of sodium atoms (representing the macroscopic Schrodingcr-cat) falling freely 
under gravity are plotted for varying wave packet separation d. In each curve, time (in 
sec) is plotted along horizontal direction and the coherent TOF distribution H(t) of 
BEC Schrodinger-cat is plotted (in sec -1 ) along vertical axis. The detector is located 
at a distance z = H = — 1 cm and (Tq = 1 /im. 



parameter values. We can see from the oscillatory factor 5 (1221) in (jTTI) that the number 
of oscillations and hence the number of fringe increases in the TOF distribution IT(t) 
( |23|) as one increase the separation d. The interference effect arises mainly because 
of two factors: one is the temporal overlap P 12 (z = H,t) J2DH and the other is the 
oscillatory factor 5. When d is very small, the overlap P\i{z = H,t) is very high, but 
the oscillatory factor 5 becomes small. As a consequence, the oscillation frequency is 
too slow or the oscillation period is too large, and we do not see any oscillatory effect 
in the temporal overlap region of the wave packets. Number of oscillations increases as 
one increases d, but again after a certain value (d > 400/im) of separation there will be 
no interference as the overlap P\i{z = H,t) becomes very small in that case. 
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Figure 3. (Color online) Quantum TOF distributions H(t) for the coherently splitted 
BEC of sodium atoms (representing the macroscopic Schrodinger-cat) falling freely 
under gravity are plotted for varying atomic masses. In each curve time (in sec) is 
plotted along horizontal direction and the TOF signal H(t) is plotted (in sec -1 ) along 
vertical axis. The detector is located at z = H = — 1 cm with d — 50 /im and 
(To = 1 M m - 



From Fig. [21 we see that the interference pattern in quantum TOF signal gradually 
disappears as one increases the mass m of the atoms. Figure H] shows the quantum TOF 
distribution for different values of wave packet width from cro = 1 /im to cr = 6 fim. It 
is clear from FigJH that the number of fringe and the contrast of interference pattern in 
quantum TOF distribution decreases as one increases the value of the initial widths (cxo) 
of the wave packets. Nevertheless, it is possible to see the interference for a larger value 
of o"o. For example, if one chooses <To = 10 /im, then to observe good interference pattern 
(with good contrast and having considerable number of fringe) in TI(t) , the separation 
d needs to be considered in the range of 50 /im to 250 /im, with the detector placed at 
a longer distance (H = — 100 cm) for a fixed mass of sodium atoms. 

The interference in H(t) is sensitive to the parameters o"o and the atomic mass m, the 
detector location H and the separation d. We repeat here that the interference in II(t) 
arises mainly because of the temporal overlap P 12 (z = H,t) flUD and the oscillatory 
factor 5 (|22|) . To increase the temporal overlap Pi2(z = H,t), one has to find the 
condition under which the spreading of the wave packet increases: small <To, lighter 
mass atoms, distant detector location (large H) will be helpful in this regard to enhance 
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Figure 4. (Color online) Quantum TOF distributions H(t) for the coherently splitted 
BEC of sodium atoms (representing the macroscopic Schrodingcr-cat) falling freely 
under gravity are plotted for varying wave packet width cto- In each curve, time (in 
sec) is plotted along horizontal direction and the coherent TOF distribution n(t) of 
BEC Schrbdinger-cat is plotted (in sec -1 ) along vertical axis. The detector is located 
at a distance z = H = — 1 cm and d = 50 fim. 



this effect. The oscillatory factor 5 can be increased either by reducing the value of do, 
or by increasing the parameters d and H. Actually, when one considers higher values 
of the parameter o"o, then the temporal overlap P\2{z = H,t) and the oscillatory factor 
S both decrease. This is because for larger values of the parameters <7o ( or mass m), 
the spreading effect ( 12TT) and hence the temporal overlap Pu(z = H,t) becomes small. 
As a result, the wave packets try to localize (in time as well as in space) more strongly 
causing the interference effect to be small. Also, for higher values of <To, the oscillatory 
factor 5 will be too small due to the presence of <Tq in the denominator of 5 (|22|) . Then 
one has to allow the BEC to travel a longer distance (by increasing H) to develope 
some temporal overlap of the wave packets, and also increasing H helps us to increase S 
([22]) . For higher values of cr , the parameter S should also be increased by increasing the 
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Figure 5. (Color online) Quantum TOF distribution FE(t) for the coherently splitted 
BEC of sodium atoms (representing the macroscopic Schrodinger-cat) is plotted in the 
absence of gravity (g = 0). Time (in sec) is plotted along horizontal direction and the 
coherent TOF distribution H(t) of BEC Schrodinger-cat is plotted (in sec -1 ) along 
vertical axis. The detector is located at a distance z = H = — 1 cm with ao = 1 /«m 
and d = 20 //m. 



value of the separation d, keeping in mind that there is a considerable temporal overlap 
Pu(z = H,t). The temporal overlap gets reduced if one increases the separation d 
too much. So, even if there is a delicate choice of the parameters, one can observe the 
interference for a wide range of parameter values. 

It is significant to mention here that gravity plays an important role in our setup to 
observe the interference in the TOF distribution. In Figl5]we plot the time distribution 
n(t) for g = (no gravity) and compare it with Fig. where we plot IT(t) in the 
presence of gravity with the parameter values same as that of Figj5l We see that the 
magnitude of U.(t) (in the absence of gravity) is roughly 10 5 times smaller than that 
obtained for gravitational free fall case for z — H — — 1 cm. The reason for this is that 
the magnitudes of J\{z = H,t) and Ji{z = H,t) become very small (roughly 10 5 times) 
in the absence of gravity. This magnitude becomes 10 6 times smaller if we consider the 
detector location at z = —10 cm. Actually, in the absence of gravity, there will be free 
particle motion and expansion of the wave packets in every direction. So, if one tries 
to observe the interference in the quantum TOF distribution in the absence of gravity, 
the intensity of that interference pattern will be too faint to be observed as only a small 
fraction of the condensate atoms will arrive at the detector. Hence, in our setup, gravity 
plays an important role which helps to pull down the condensate towards the detection 
plane. 
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4. Summary and Conclusion 

To summarize, in this work we propose a scheme to experimentally observe matter- 
wave interference in the time domain, specifically in the TOF (arrival-time) distribution 
using atomic BEC. This experimentally testable scheme has the potential to empirically 
resolve ambiguities inherent in the theoretical formulations of the quantum arrival 
time distribution. Here we use the probability current density approach to calculate 
the quantum TOF distributions for atomic BEC Schrodinger cat represented by 
superposition of macroscopically separated wave packets in space. Our definition of the 
quantum TOF distribution in terms of the modulus of the probability current density is 
particularly motivated from the equation of continuity, and other physical considerations 

discussed in the literature [m[IEl[T3l2dl2DI23l2SI2SI25]. This a PP roach also provides 
a proper classical limit, as the interference and hence the coherence in the quantum TOF 
signal disappears in the large- mass limit. We repeat that there is no classical analogue 
of this TOF distribution H(t) and this is purely a quantum distribution where we see 
the matter-wave interference in the quantum TOF signal. Hence, it will be interesting 
to see if our prediction of interference in time domain (TOF distribution) can be verified 
in actual experiments using modern interferometry techniques and sophisticated TOF 
methods. 
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